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Abstract 

Existence check of non-trivial, stationary axisymmetric black hole solutions 

in Brans-Dicke theory of gravity in different direction from those of Penrose, 

Thorne and Dykla, and Hawking is performed. Namely, working directly 

with the known explicit spacetime solutions in Brans-Dicke theory, it is found 

that non-trivial Kerr-Newman-type black hole solutions different from general 

relativistic solutions could occur for the generic Brans-Dicke parameter values 

— 5/2 < lo < —3/2. Finally, issues like whether these new black holes carry 

scalar hair and can really arise in nature and if they can, what the associated 

physical implications would be are discussed carefully. 
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I. Introduction 

Of all the alternative theories of classical gravity to Einstein's general relativity, perhaps 
the Brans-Dicke (BD) theory [1] is the most studied and hence the best-known. This theory 
can be thought of as a minimal extension of general relativity designed to properly accomo- 
date both Mach's principle [2] and Dirac's large number hypothesis [2]. Namely, the theory 
employs the viewpoint in which the Newton's constant G is allowed to vary with space and 
time and can be written in terms of a scalar ("BD scalar") field as G — 1/<E>. In this work, 
we are interested in the existence of exact solutions to the BD field equations that can 
describe rotating, charged black hole spacetimes and their detailed structure. And if there 
are, we would like to know whether they are non-trivial ones different from general relativis- 
ts black hole solutions. As is well-known, even in Einstein's general relativity, to find the 
exact solutions to the highly non-linear Einstein field equations is a formidable task. For 
this reason, algorithms generating exact, new solutions from the known solutions of simpler 
situations have been actively looked for and actually quite a few were found. In BD theory 
of gravity, the field equations are even more complex and thus it is natural to seek similar 
algorithms generating exact solutions from the already known simpler solutions either of 
the BD theory or of the conventional Einstein gravity. To the best of our knowledge, meth- 
ods thus far discovered along this line includes those of Janis et al., Buchdahl, Mcintosh, 
Tupper, Tiwari and Nayak, and Singh and Rai [3]. In particular, Tiwari and Nayak [3] pro- 
posed an algorithm that allows us to generate stationary, axisymmetric solutions in vacuum 
BD theory from the known Kerr solution [6] in vacuum Einstein theory and later on Singh 
and Rai [3] generalized this method to the one that generates stationary, axisymmetric, 
charged solutions in BD-Maxwell theory from the known Kerr-Newman (KN) solution [6] in 
Einstein-Maxwell theory. Thus in the present work, we shall take, as the Kerr-Newman-type 
solutions in BD-Maxwell theory (henceforth "BDKN" solutions), the ones constructed by 
Singh and Rai to explore if it can descibe non-trivial black hole spacetimes different from 
those described by the standard KN solution in Einstein-Maxwell theory. 
II. Non-trivial BDKN black hole solutions 
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We begin by briefly reviewing the algorithm proposed first by Tiwari and Nayak and gener- 
alized later by Singh and Rai. Consider the BD-Maxwell theory described by the action 
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where <£> is the BD scalar field and uj is the generic parameter of the theory. Extremizing this 
action then with respect to the metric g^, the BD scalar field $, and the Maxwell gauge 
field (with the field strength F^ v = V^A U — V U A^) yields the classical field equations 
given respectively by 

(2) 

7^V U <5> - \g^V a ® V a $ ) + ^(V„V„$ - <^V Q V Q $ ) 
and 

V Q V a $ = Q , Tf A = 0, V,F^ = 0, V,F^ = 

with the last equation being the Bianchi identity and F^ u = |e^F a/3 . And the Einstein- 
Maxwell theory is the uj — > oo limit of this BD-Maxwell theory Note that in the action 
and hence in the classical field equations, there are no direct interactions between the BD 
scalar field $ and the ordinary matter, i.e., the Maxwell gauge field A^. Indeed this is the 
essential feature of the BD scalar field $ that distinguishes it from "dilaton" fields in other 
scalar-tensor theories such as Kaluza-Klein theories or low-energy effective string theories 
where the dilaton-matter couplings generically occur as a result of dimensional reduction. 
(Here we would like to stress that we shall work in the context of original BD theory format 
not some conformal transformation of it.) As a matter of fact, it is the original spirit [1] of 
BD theory of gravity in which the BD scalar field <E> is prescribed to remain strictly massless 
by forbidding its direct interaction with matter fields. Now the algorithm of Tiwari and 
Nayak, and Singh and Rai goes as follows. Let the metric for a stationary, axisymmetric, 
charged solution to Einstein-Maxwell field equations take the form 



d S 2 = _ e W E ( dt + WEd( pf + e 2(k E -U E )^ dx lj2 + (^2)2] + ^2 e 
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while the metric for a stationary, axisymmetric, charged solution to BD-Maxwell field equa- 
tions be 



where [/, W, fc and /i are functions of re 1 and x 2 only. The significance of the choice of the 
metric in this form has been thoroughly discussed by Matzner and Misner [4] and Misra and 
Pandey [5] . Tiwari and Nayak, and Singh and Rai first wrote down the Einstein-Maxwell and 
BD-Maxwell field equations for the choice of metrics in eq.(3) and (4) respectively. Com- 
paring the two sets of field equations closely, they realized that stationary, axisymmetric 
solutions of the BD-Maxwell field equations are obtainable from those of Einstein-Maxwell 
field equations provided certain relations between metric functions hold. 
That is, if (We, &e, Ue, He, A?) form a stationary, axisymmetric solution to the Einstein- 
Maxwell field equations for the metric in eq.(3), then a corresponding stationary, axisym- 
metric solution to the BD-Maxwell field equations for the metric in eq.(4) is given by 
(W B d, k BD , U BD , h BD , A^ D ) where 



Now what remains is to apply this method to obtain the Kerr- Newman-type solution in 
BD-Maxwell theory (BDKN solution) from the known Kerr-Newman (KN) solution [6] in 
Einstein-Maxwell theory. And to do so, one needs some preparation which involves casting 
the KN solution given in Boyer-Lindquist coordinates [7] (t, r, 6, 0) in the metric form in 
eq.(3) by performing a coordinate transformation (of r alone) suggested by Misra and Pandey 
[5]. Namely, we start with the KN solution written in Boyer-Lindquist coordinates 



W B d = W E , k BD = k E , U B d = U E - - log $, 
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where E = r 2 + a 2 cos 2 and A = r 2 — 2Mr + a 2 + e 2 with M, a, and e denoting the ADM 
mass, angular momentum per unit mass, and the electric charge respectively. Consider now 
the transformation of the radial coordinate introduced by Misra and Pandey [5] 



o t , (M 2 -a 2 -e 2 

r = e R + M + - 
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which gives dr 2 /A = dR 2 . Then the KN solution can now be cast in the form in eq.(3), i.e., 

j9 r A — a 2 sin 2 0, , , TTT , _ r ,„ 9 , A sin 2 6> , , 9l . . 

ds 2 = - [ (dt + Wd(f)) 2 + £ [d# 2 + rfi? 2 + ——rf^ 2 (8) 

L A — a z sin 

with now £ = L 2 + a 2 cos 2 and A = L 2 — 2ML + a 2 + e 2 where we set, as a short-hand 

notation, L = e R + M + ^ M —. ~ e ^ e~ R . Now we can read off the metric components as 



2Ue _ A - a 2 sin 2 _ a sin 2 0(L 2 + a 2 - A) 

6 B_i £ J ' £ ~ A - a 2 sin 2 ' 
e 2k E = ( A _ a 2 sin 2 g ^ h 2 E = A sin 2 0. 
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Then using the rule in eq. (5) in the algorithm by Tiwari and Nayak, and Singh and Rai, we 
can now construct BDKN solution in BD-Maxwell theory as 
'L 2 + a 2 cos 2 - 2ML + e 2 * 
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+ ( L 2 + a 2 _ 2 ML + e 2)(2-+D/(2- + 3) sin 2(2a, + l)/(2a, + 3) g ^ 

0) = (L 2 + a 2 - 2 ML + e 2 ) 2 ^ +3 ) sin 4 /^+ 3 ) 0, 

A ^- L 2 + lLs 2 ^- a ^ 9 ^ 
Then by transforming back to the standard Boyer-Lindquist coordinates using eq.(7), we 
finally arrive at the BDKN solution in Boyer-Lindquist coordinates given by 

ds 2 = A- 2 ^ 2 ^ sin- 4 /^+ 3 ) \- r— - | dt- - z__^__ ^L ( j i(Jo 
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Note that as uj — > oo, this BDKN solution goes over to the standard KN solution as it should 
since the uj — > oo limit of BD theory is the Einstein gravity. Now that we have an exact, 
electrovac, stationary axisymmetric solution to the BD-Maxwell theory. Then it is natural to 
ask if this BDKN solution can describe a black hole spacetime resulted from a gravitational 
collapse. And if it can, furthermore one might be curious whether this BDKN solution or its 
three special cases (i.e., BD-Schwarzschild, BD-Reissner-Nordstrom or BD-Kerr solutions) 
could describe non-trivial black hole spacetimes which are different from those described by 
their general relativistic counterparts. Indeed, questions of this sort had been raised long 
ago, and actually Penrose [8] conjectured that even in BD theory of gravity, the relativistic 
gravitational collapse in three spatial dimensions would produce black holes identical to those 
in general relativity. And this conjecture received some support from the work of Thorne 
and Dykla [8] in which they presented four pieces of evidence in favor of the conjecture 
by employing mainly the "large-u;" expansion scheme (recall that in the limit uj — > oo, the 
BD theory goes over to the general relativity). As Thorne and Dykla mentioned in their 
work, however, the conjecture of Penrose was not fully proved since detailed analysis of the 
collapse with arbitrary, finite values of the generic BD parameter uj is needed. In this regard, 
we now seem to be in a better shape toward the serious investigation on the validity of the 
conjecture since we have an exact, stationary axisymmetric solution possessing arbitrary 
uj values which was not available at the time. Therefore we begin with the bottomline 
qualification for the BDKN solution in BD-Maxwell theory to describe a rotating, charged 
black hole spacetime, namely the possible occurrence of non-singular event horizon. Then 
along this line, perhaps the most natural first step is to ask under what circumstances the 
Killing horizons develop. Just like the KN solution in Einstein-Maxwell theory, this BDKN 
solution is stationary and axisymmetric and hence possesses the time translational Killing 
field £ M = (d/dty and the rotational Killing field ^ = (d/d4>Y correspondingly and it is 
their linear combination, — £ M +fiff"0^ which is normal to the Killing horizons, if any (here 
Qh denotes the angular velocity of the Killing horizon). And if Killing horizons are present, 
they occur at points where x^ becomes null which turn out to be zeroes of /\( 2u} + 1 )/( 2uJ + 3 ) — o. 

6 



Thus first for uo = —1/2, obviously no horizon occurs. Next for (§^§) > 0, i.e., for 



uj < —3/2 or u) > —1/2, two Killing horizons occur at r± = M± (M 2 



3 2U/2 



(provided 



M 2 > a 2 + e 2 ) which are precisely the same locations as those of Killing horizons of KN 
black holes in Einstein-Maxwell theory. Now since the formation of horizons appears to be 
possible, next we investigate their nature. And to this end, we examine behaviors of the 
invariant curvature polynomials such as R, R^ V R^ U and R^apR^ 11 "* 13 , the surface gravity 
k and the energy density of the BD scalar field T^f^^ on these candidates for Killing 
horizons. And as we mentioned above, since the bottomline qualification for the black hole 
interpretation of BDKN solution is the regularity of the horizon candidate, we begin with 
the examination of behavior of invariant curvature polynomials on the horizon candidate at 
which A = 0. First, the curvature scalar is calculated to be 
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As was the case with KN black hole solutions, it also blows up at £ = (i.e., r = 0, = tt/2) 
indicating that the BDKN black hole solution also has the curvature singularity with the 
same "ring" structure. The direct computation of the other two curvature polynomials, 
i.e., the Ricci square R^ U R^ U and the Kretschmann curvature invariant R^uapR^^ for this 
BDKN solution is a formidable job. But a close inspection reveals that indeed we can save 
considerable amount of labor. Namely, consider now the Brans-Dicke-Schwarzschild (BDS) 
spacetime solution that can be obtained by setting a = e = in the BDKN solution in 
eq.(ll) 

2AT 



ds 2 = A- 2 ^ + V sin-V^+s) e 

+ A 2/(2^+3) sin 4/(2^+3) g 
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$ (r,0) = A 2 / (2 ^ sin 4 /^ +3 )# 



where now A = r(r — 2M). A remarkable feature of this BDS solution is the fact that, unlike 
the Schwarzschild solution in general relativity, the spacetime it describes is static (i.e., non- 



rotating) but not spherically-symmetric. Thus computing R^W and R^apR^ 013 for this 
BDS solution and examining their behaviors on the horizon candidate at which A = would 
be sufficient to envisage the possibility of regular event horizon for both BDS and BDKN 
solutions. In addition, another noticeable characteristic of both BDKN and BDS spacetime 
solutions is that they have possible coordinate singularities not only at the outer event 
horizon where A = but also along the symmetry axis 9 = 0, ir. Thus in order to explore 
the nature of this metric singularity along the symmetry axis, the computation of invariant 
curvature polynomials looks necessary. The result of the computation of the Ricci square and 
the Kretschmann curvature invariant for this BDS solution is given in the appendix. And it 
is a straightforward matter to realize that the two invariant curvature polynomials R^ U R^ V 
and R^apR^ 013 given in the appendix and the curvature scalar R given above in eq.(12) 
become finite both on the horizon candidate at which A = and along the symmetry axis 
9 = 0, 7r provided the generic BD ^-parameter takes values in the range —5/2 < uo < —3/2. 
Next, we turn to the computation of surface gravity at the Killing horizons, k±. Generally, 
the surface gravity k is defined in a gravity theory- independent manner as follows. Since 
the horizon is a null surface, there we have x^Xn — where x^ is the Killing field normal 
to the horizon we defined above. This implies that V^ix^Xv) is also normal to the horizon. 
Thus on the horizon, there exists a function k such that 

VWXu) = -2^X M or x V V„X» = 
from which it can be derived that 

^ = ^(VY)M (14) 

where the evaluation on the horizon is understood. Now, for the non-trivial BDKN black 
hole solution at hand, a straightforward albeit somewhat tedious calculation yields 

Lastly, the energy density of the BD scalar field is computed using eq.(2) as 
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It is interesting to note that this energy density of the BD scalar field also blows up at the 
curvature singularity £ = 0. 

Now (i) for lo -> oo, R = 0, re± = (r± - r T )/2(r| + a 2 ) and T^ff = on the 
surfaces r = r±. This is an anticipated result since this is the correct KN black hole 
limit in Einstein-Maxwell theory, (ii) Next for oo > u > —1/2, on the surfaces r = r±, 
(R, R^Ri", R^apR^) -> oo, k± -> oo and T^ D ^ U -> oo with $(r±,0) = 0. This 
indicates that the surfaces r = r ± are singular and fail to act as horizons and hence the 
corresponding metric cannot describe a black hole spacetime. (iii) Finally for oo < —3/2, or 
more precisely for —5/2 < uo < —3/2, on the surfaces r = r±, (.R, R^R^ , R ilva pR iiVali ) = 
(or const, particularly for a; = -5/2), « ± = and T^ V D ^ V = with $(r±, 0) — > oo. Namely 
the curvature invariants are finite, surface gravity is zero and the BD scalar field satisfies the 
weak energy condition although its value diverges there. (Here, infinite value of $ indicates 
that the effective Newton's constant tends to zero.) Thus in this range of the u- values, the 
surfaces r = r± may act as regular Killing horizons and hence the corresponding BDKN 
metric solution appears to describe non-trivial black hole spacetimes different from those in 
Einstein-Maxwell theory. In particular for cu = —5/2, the corresponding non-trivial BDKN 
black hole solution singles out with a relatively simple form given by 
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(17) 
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Lastly, for the rest of the upvalues, i.e., for lo < —5/2, the surfaces r = r± are singular 



horizons on which (R, R^R^, R^ a pR^ aP ) -> oo, k + = and T^ D ^ U = and for 
—3/2 < < —1/2, no Killing horizon develops and thus the curvature singularity at S = 
is naked. Therefore it now appears that for the values of the generic BD cj-parameter in the 
limited range —5/2 < uo < —3/2, the BDKN solution in BD-Maxwell theory may describe 
non-trivial black hole spacetimes. 
III. Nature of BDKN black holes 

Now that we have non-trivial BDKN black hole solutions. It seems then natural to explore its 
thermodynamics and causal structure in some more detail. Firstly, these BDKN black hole 
solutions have vanishing surface gravity at the event horizon, k+ — and hence zero Hawking 
temperature, T H = k + /2ti = 0. In other words, they do not radiate and hence are completely 
"dark and cold". Certainly, this is a very bizzare feature in sharp contrast to evaporating 
black holes in general relativity. Next, we turn to their causal structure. As noted earlier, 
the two Killing horizons, i.e., the outer event horizon and the inner Cauchy horizon turn 
out to occur precisely at the same locations (i.e., same coordinate distances) as those of KN 
black hole solutions in Einstein-Maxwell theory, i.e., at r± — M ± (M 2 — a 2 — e 2 ) 1 / 2 . Also 
it is interesting to note that the proper area of the event horizon at r = r + , 

A = [ dOd^geeg^) 112 = 47r(r 2 + a 2 ) (18) 

is again exactly the same as that of standard KN black hole spacetime. In addition, its 
angular velocity at the event horizon coincides with that of standard KN solution as well 

- WEb(r + ) = = -W E \r + ). (19) 

Next, observe that the norm of the time translational Killing field 



^ = gtt = _ A -2/(2u,+3) sin -4/(2u,+3) q 



A - a 2 sin 2 9 



(20) 



goes like negative (r_ < r < r+) — > positive (r+ < r < r s ) — > negative (r > r s ) with 
r s = M + (M 2 — a 2 cos 2 # — e 2 ) 1 / 2 > r + being the larger root of indicating that ^ 

behaves as timelike — > spacelike — > timelike correspondingly. And particularly the region in 
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which £ M stays spacelike extends outside hole's event horizon. This region is the so-called 
"ergoregion" and its outer boundary on which £ M becomes null, i.e., r = r s is called "static 
limit" since inside of which no observer can possibly remain static. Thus if we recall the 
location of the static limit in standard KN black hole solution, we can realize that even 
the locations of ergoregions in two black hole spacetimes, KN and BDKN, are the same as 
well. Namely in two theories, i.e., the BD-Maxwell theory and the Einstein-Maxwell theory, 
rotating, charged black hole solutions turn out to possess identical causal structure (i.e., the 
locations of ring singularities, two Killing horizons and static limits are the same) and hence 
exhibit the same global topology. Thus actually what distinguishes the BDKN black hole 
spacetime from its general relativity's counterpart, i.e., the KN black hole is the local geom- 
etry alone such as the curvature characterized by the specific ui- values, —5/2 < uj < —3/2. 
At this point, perhaps it is relevant to mention the behavior of the BD scalar field which 
plays unique role only in BD theory of gravity. Independently of Penrose [8] and of Thorne 
and Dykla [8], Hawking [9] also explored the possible existence of black hole solutions in BD 
theory and put forward a theorem which states that stationary black holes in BD theory 
are identical to those in general relativity. To be a little more concrete, Hawking extended 
some of his theorems for general relativistic black holes to BD theory and showed that any 
object collapsing to a black hole in BD theory must settle into final equilibrium state which 
is either Schwarzschild or Kerr spacetime. And in doing so, he "assumed" that the BD scalar 
field $ satisfies the weak energy condition and is constant outside the black hole. Therefore 
now one may be puzzled as we realized in the present work that non-trivial BDKN black 
hole solutions different from general relativistic KN solution could exist in seemingly con- 
tradiction to Hawking's theorem. There is, however, no contradiction. Hawking deduced 
the theorem by manipulating the BD field equations and most crucially "assuming" the 
strict conditions on the BD scalar field stated above but not by working with an explicit 
spacetime solution which was not available at the time. In the present work, however, we 
investigated closely the known, explicit stationary axisymmetric solution in BD theory. And 
in particular, when the BDKN solutions can describe non-trivial black hole spacetimes for 
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the BD oparametr values —5/2 < cu < —3/2, the accompanying BD scalar field solution 
$(r, 9) = A 2/(2w+3) s in 4/(2w+3) 9 turns out not to be a constant field outside the event horizon 
at r = r+. Besides, the energy density of this BD scalar field, T^J } ^ ti ^ v whose explicit 
form was given earlier in eq.(16) does not strictly obey the weak energy condition for all r. 
Namely the value of T^ U D ^^ does not remain no n- negative for all r. Rather, its value and 
hence the signature changes from point to point. In short, the Hawking's theorem simply 
cannot be applied to the present situation and hence the results of the present study needs 
not be restricted by Hawking's theorem. At this point, it seems appropriate to ask whether 
the non-trivial BDKN black hole solution studied in the present work can be viewed as a 
counterexample to the no-hair theorem of black holes. In the loose sense, one may think of 
the non-trivial behavior of the BD scalar field outside the event horizon as indicating the 
appearance of "scalar hair". Here, however, we need to be more precise with the nature 
of no-hair theorem. Following Bizon [11], for instance, a certain theory is said to allow a 
hairy black hole solution if there is a need to specify quantities other than conserved charges 
defined at asymptotic infinity such as the mass, angular momentum and the electric charge 
in order to characterize comletely a stationary black hole solution within that theory. Thus 
in this stricter sense, the non-trivial BDKN black hole solution studied in the present work 
does not constitute a hairy black hole solution since both the metric and BD scalar field 
solutions in eq.(ll) are specified completely by the ADM mass M, angular momentum per 
unit mass a, and the electric charge e only and no other quantities. 
IV. Discussions 

Before we address the physical implication of the non-trivial BDKN black hole solution 
found in the present work, we would like to comment on a technical issue, i.e., the diver- 
gent behavior of the BD scalar field solution on the horizon. And in relation to this, it is 
interesting to note that our BDKN black hole solution shares two peculiar features, i.e., the 
divergent behavior of the scalar field on the horizon and the null Hawking radiation, with 
the well-known Bekenstein black hole solution in Einstein-conformal scalar field theory [12]. 
Namely, using a suitable solution generation technique, long ago, Bekenstein constructed 
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a static, spherically-symmetric black hole solution in which the metric part corresponds 
to that of extreme Reissner-Nordstrom (RN) black hole and thus represents non-radiating 
black hole spacetime and the conformal scalar field solution diverges on the horizon. There- 
fore for direct and parallel comparison between the two solutions, it seems appropriate to 
take the BDS solution analyzed in detail in the appendix. The two theories, of course, 
have completely different nature and motivations. The Einstein-conformal scalar field the- 
ory has been devised guided mainly by a particular (Weyl) symmetry and the scalar field 
there is supposed to describe a matter. The BD theory, on the other hand, is an alterna- 
tive theory to Einstein gravity and the BD scalar field here represents a spacetime-varying 
effective Newton's constant, not a matter. Thus the divergent behavior of the scalar field 
in Einstein-conformal scalar field theory could be disastrous but that of the BD scalar field 
in BD theory essentially represents the vanishing effective Newton's constant in a certain 
region of spacetime. Besides, since the energy density of the BD scalar field T^J 3 ^^ given 
in eq.(16) vanishes and hence satisfies the weak energy condition on the horizon at which 
A = (of course for —5/2 < uj < —3/2), we do not worry too much about the divergent 
behavior of the BD scalar field there. Another interesting contrast is that the metric solution 
of Bekenstein black hole spacetime there corresponds to a familiar extreme RN metric which 
is static and spherically-symmetric whereas the BDS metric solution here exhibits a remark- 
able feature that it is static (i.e., non-rotating) but not spherically-symmetric as we pointed 
out earlier. Besides, the Bekenstein solution represents a "hairy" black hole [13] whereas our 
BDS black hole solution carries no hair as we mentioned above. Now, the point of central 
interest we would like to make is about the issue raised recently by Sudarsky and Zannias 
[13]. To be a little more concrete, they showed that the divergent behavior of the conformal 
scalar field solution on the horizon essentially leads the Bekenstein black hole solution to fail 
to satisfy Einstein field equations particularly on the horizon. They, thus, concluded that 
the Bekenstein solution cannot be considered as a genuine black hole solution and therefore 
the black hole no-hair theorem is saved. And as a manifest evidence for their argument 
against the black hole interpretation, Sudarsky and Zannias demonstrated that by working 
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in Eddington-Finkelstein null coordinates, parts of the Einstein field equations can be shown 
not to hold in a rigorous sense as the left-hand side of the equation, say R^ v vanishes on 
the horizon (since the metric is that of extreme RN) while the right-hand side, namely the 
energy-momentum tensor of the conformal scalar field, 87r[T Mi , — is ill-defined on 

the horizon as the conformal scalar field diverges there. Thus with this in mind, now we 
consider the validity of our BDS solution on the horizon. As one can see in the appendix, the 
5 non-vanishing components of BD field equations, i.e., tt, rr, r9, 99 and 00 parts appear 
to hold perfectly. In particular, on the horizon and for —5/2 < uj < —3/2, tt and 00 parts 
hold as "0 = 0". The other 3 parts, however, hold as "oo = oo". Namely, in these 3 equa- 
tions, not only the energy-momentum tensor of the BD scalar field on the right-hand side 
of BD field equations but also the Ricci tensor components on the left-hand side diverges 
in exactly the same manner on the horizon. Indeed precisely these 3 equations are the ones 
we need to be careful with. Here, however, we must say that a naive attempt toward the 
validity check of the BDS solution on the horizon in exactly the same way as Sudarsky and 
Zannias did for Bekenstein solution seems to be obscured. In fact, the demonstration of Su- 
darsky and Zannias was successful largely because the simple transformation from spherical 
to Eddington-Finkelstein null coordinates for the Bekenstein solution was available. Indeed, 
the virtue of null coordinates is that in terms of which g rr = and hence / M = (d/dr)^ 
becomes a smooth null vector field such that the quantity R^l^V can be shown to be finite 
(zero) on the horizon while the right-hand side of Einstein equation, 87r[T^ — g ^JT^ / 2]IH V 
is ill-defined. In contrast, however, taking a Eddington-Finkelstein-type null coordinates for 
the BDS solution is not so obvious as in the case of static, spherically-symmetric black hole 
solutions in Einstein theory and is indeed practically awkward. It is essentially due to the 
peculiar feature of BDS metric solution which is static but not spherically-symmetric as has 
been stressed earlier. Thus in the present work, we do not pursue validity check of the BDS 
solution on the horizon in this direction and leave it as an issue for future investigation. As 
a result, for the 3 parts of the BD field equations yielding "oo = oo", no definite statement 
can be made yet concerning whether or not they are valid, i.e., these 3 equations are really 
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satisfied on the horizon. Nevertheless, one interesting point we can make is that if the BDS 
solution does represent a genuine black hole spacetime, then the black hole no-hair theorem 
appears to survive even in the BD theory of gravity since the BDS (and BDKN as well) 
spacetime is not a hairy black hole solution as we discussed earlier. 

It seems that now the most relevant question to ask is ; if they are genuine, would these 
non-trivial BDKN black hole spacetimes in BD theory of gravity exhibiting bizzare features 
such as null radiation really arise in nature? Of course this question needs to be answered 
very carefully and honestly and the answer for now does not seem to be in the affirmative. 
Firstly, from field theory's viewpoint, the generic BD theory ^-parameter has to be "posi- 
tive" in order for the BD scalar field <3> to have canonical (positive-definite) kinetic energy 
as can be seen in the BD gravity theory action given in eq.(l). Secondly, it is well-known 
that the BD gravity theory is in reasonable accord with all available observations and ex- 
periments thus far provided \uj\ > 500 [10]. Since both these constraints on the values of 
the cu-parameter seem to rule out the range — 5/2 < uo < — 3/2 in which the BDKN solution 
could describe non-trivial black hole spacetimes, for now it seems fair to say that these non- 
trivial BDKN black hole spacetimes different from their general relativistic counterparts are 
unlikely to arise in nature. This, however, may not be the end of the story. As we have seen 
in this work, the energy density of the explicit BD scalar field solution (which essentially 
consists of its kinetic energy) turns out not to satisfy the weak energy condition irrespective 
of the us- value. Perhaps this implies that we may abandon the "canonical kinetic energy" 
condition on the BD scalar field and allow negative-^ values. Moreover, the lower bound 
> 500 may be relaxed considerably with the advances in technology associated with 
astronomical observations and astrophysical experiments. Thus perhaps it might be wise 
to keep the possibility of non-trivial BDKN black holes alive. As a matter of fact, there 
is another type of possibility of greater physical significance and relevance. Note that the 
generic BD cu-parameter is a kind of coupling constant appearing in the BD gravity action. 
Thus in principle, it should be considered as a "running" coupling constant as a result of 
renormalization in the quantum gravity context. And its scale-dependent behavior can be 
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envisaged as follows. In the BD gravity action given in eq.(l), the term ~ tu(V a &V a $/$), 
like other terms in the action, should be finite. Thus large-w indicates the regime where 
the BD scalar field $ remains nearly constant which corresponds to the large-scale present 
universe limit (in which the BD theory goes over to the general relativity). On the other 
hand, small-w indicates the regime where the BD scalar field varies sizably with space and 
time which would presumably correspond to the small-scale early universe limit. Thus if 
we are willing to accept the BD theory as a "better" effective theory of quantum gravity 
than general relativity to describe the entire stages (scales) of the universe evolution, then 
at early times when the value of uj was small such as —5/2 < uj < —3/2, the non-trivial 
BDKN black holes like the ones studied in this work would have had a chance to form. 
These "primordial" black holes, unlike their general relativistic counterparts, however, do 
not evaporate as we discussed earlier. Thus it can be speculated that they might still hide 
somewhere in the dark side of the space today as a possible constituent of the cold dark 
matter. After all, new discoveries can be made when we keep our minds as well as eyes open. 



Acknowledgements 

This work was supported by grant of Post-doc. Program at Kyungpook National Uni- 
versity (1998). 



Appendix : Computation of invariant curvature polynomials 

In this appendix, we shall explicitly write down the BD field equations satisfied by the 
BDS solution given in eq.(13) and then provide the result of the computation of its invariant 
curvature polynomials such as the Ricci square R^ U R^ V and the Kretschmann curvature 
invariant R^R^ . 
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First we start with the BD field equations. The tensor part of the vacuum BD field equations 
that this BDS solution satisfies, 



OJ 



1 



R — v <f>V $ -I V V $ 



(21) 



has the following 5 non-vanishing components. Here we show that actually they are all 
satisfied by the BDS solution in the sense that the explicit computation of the left-hand 
side, i.e., and the right-hand side, i.e., V^V^ + ^V M Vj,<& performed with the BDS 
solution given in eq.(13) precisely agree. 



Rtt = |^V t $V t $ + -j-ViV t $ 



{2oj + 3) 



(2w + 3) 



(r - Mf - M(r - M) + 



(2oj + 3) 



cot 2 OA 



oj 1 
R rr = _v r $V r $ + -v r v r $ 



1 



A 2 (2w + 3) 



(2oj + 3) 



(r - M) 2 + M(r — M) + < 1 + 



(2w + 3) 



cot 2 \ A 



oj 1 



1 n 4 
— COt 6*7- r 

A (2w + 3) 



4w 



(2w + 3) 



(r — M) — (r — 2M) 



(22) 



oj 1 



A (2w + 3) 



(2w + 3) 
1 



(r - Mf + (r - M)(r - 2M) - <^ 1 



2oj - 1 
277+3 



cot 2 y a 



ZlA- 4 /(2^+3) sin -8/(2o;+3) g 

A 



{ (2^T3) (r - M)2 - (r - M){T - 2M) ] ^ " - ~" " A 



(2w + 3) Ll(2w + 3) v v ~~ /v "') V2^ + 3, 

And the scalar part of the vacuum BD field equations can easily be seen to hold. 
Next we turn to the expressions for the invariant curvature polynomials. The Ricci square 
is calculated to be 



R^R^ = g^g^R^Rap 

J_ A -2(2o,+5)/(2a,+3) • -8/(2w+3) a 
r 4 



(23) 



{2oj + 3) 



(2oj + 3) 



(r - Mf - M(r - M) 



+ 



(2oj + 3) 



cot 2 6 A 
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+ 
+ 
+ 



(2a; + 3) 
4 

(2a; + 3) 
4 

(2a; + 3) 



(2a; + 3) 
2 

(2o> + 3) 
2 

(2a; + 3) 



(r - M) 2 + M(r - M) + 1 + 



cot 2 e i a 



(2a; + 3) 

(r - M) 2 + (r - M)(r - 2M) - - (fj^|) c °t 2 A 



(r-M) 2 -(r-M)(r-2M) 



/ 2a; + l \ 
1.2a; + 3/ 



cot 2 #A 



+ ± A - (2^+7)/ (2^+3) sin -2(2a,+7)/(2a,+3) g cog 2 g 



(2w + 3) 



4a; ^ 2 

(r — M) — (r — 2M) 



(2a; + 3) 



It is worth noting that this Ricci square vanishes in the limit u — > oo as it should since 
in which the BDS solution goes over to the Schwarzschild solution. And the Kretschmann 
curvature invariant is computed to be 



(24) 



_ A -2(2u,+5)/(2u,+3) sin -8/(2a,+3) q 



I { 2 {2u + 5 l (r — M) 2 — ( ^±T) M(r - M) + M 2 



(2w + 3) 2 



(2a; + 3) 



1 + 



(2a; + 3) 



cot 2 



A 



+ 4 
+ 



(2w + 3) 
2 



(r - M) + (r — 2M) 



(2a; + 3) 

+ {(£f>-" » 



(r - M ) - (r - 2M) 



(2a; + 3) 
2 



(r — M) — M 



V 2a; + 5 



1 



(2a; + 3) 



1 + 



(2a; + 3) 



(r-M)-M 



(2a; + 1) - n . 
' cot 2 6>A 



(2a; + 3) 

2 



cot 2 



A 



(2a; + 3)' 



(2a; + 3) 



(r-M)-(r- 2M) 



\2uj + 3j 



1 + 



(2a; + 3) 



cot 2 



A 



+ ^ A -2(2w+5)/(2w+3) sin -4(2w+5)/(2w+3) q 



(2w + 3) 



(r — M) + (r — 2M) 



+ 



(2a; + 3) 
2 



sin 



(2a; + 3) 



(r-M)-M 
(r-M)-(r- 2M) 



(2a; + 3) 
2 

(277+3) 



1 + 



(2w + 3) 



cos 



A 



(r-M)-M 



sin 2 0-2 



(2w + 1) 2 fl a 

cos 2 6A 



+ 2 <{ (r — M) 

'2a; + 7 
2a; + 3 
2 



+ 
+ 



(2a> + 3) 
(r-M) 



(r — M) + {r — 2M) 



sin 2 6 + 



(2w + 3) 2 
2 /2a; + 5\ 

(2a; + 3) ~ \2a; + 3/ 



sm 



A 



(2a; + 3) 
2 



(2a; + 3) 



(2a; + 3) 
(r — M) + (r — 2M) 
4 



(r-M)-(r- 2M) 
2 



2 „ , /2a; + 1\ 



sm 1 9 + 



\2a; + 3/ 



, /2a> + 5\ 2 

1 - COS^ I 

V2a; + 3/ 



A 



(2a; + 3) 



(r-M)-M 



sin 2 9 



1 + 



(2a; + 3) 



cos 2 # 



A 



(25) 



+ 2 



(2oj + 3)2 



(r - M) 2 - (r - 2Mf 



\2lo + 3J 



\2uj + 3J 



cos 



A 
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+ —&-(^+7)/^+?>) sin -2(2o;+7)/(2a;+3) COS 2 6> 

r 4 



1 



(2w + 3) 



6 V 2 
(r — M) + (r — 4M) 



(2w + 3) 



+ 



1 



(2w + 3) 



4cu 



(2w + 3) 



(r - M) + (r - 2M) 



Again, it is straightforward to check that this Kretschmann curvature invariant reduces to 
that of Schwarzschild solution, 48M 2 /r 6 in the Einstein gravity limit (i.e., as uo — > oo) as 
it should. In conclusion, a close inspection reveals that R^W 11 and R^^R^ 01 ^ are finite 
(more precisely vanishes) both on the horizon candidate at which A = and along the 
symmetry axis = 0, n provided the generic BD a;-parameter takes values in the range 
-5/2 < uj < -3/2. 
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